Topological Summation in Lattice Gauge Theory 
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Abstract. In gauge theories the field configurations often occur in distinct topological sectors. 
In a lattice regularised system with chiral fermions, these sectors can be defined by referring 
to the Atiyah-Singer Index Theorem. However, if such a model is simulated with local updates 
1 of the lattice gauge configuration, the Monte Carlo history tends to get stuck in one sector for 

^ many steps, in particular on fine lattices. Then expectation values can be measured only within 

specific sectors. Here we present a pilot study in the 2-flavour Schwinger model which explores 
methods of approximating the complete result for an observable — corresponding to a suitable 
sum over all sectors — based on numerical measurements in a few specific topological sectors. 
We also probe various procedures for an indirect evaluation of the topological susceptibility, 
^/-^ , starting from such topologically restricted measurements. 
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fNJ . 1. Topological sectors in gauge theories 



Our general framework in this article is the functional integral formulation of quantum physics 
in Euclidean space. In this setting, the set of configurations may occur in disjoint subsets, so that 
all continuously deformed configurations belong to the same subset. Such subsets are known as 
topological sectors. Continuous deformations capture all configurations in one topological sector, 
but none of any different sector (general aspects are discussed e.g. in Refs. PQ). 

The simplest example where this situation occurs is a quantum mechanical scalar particle 
moving on the circle S , with periodic boundary conditions in Euclidean time. The expectation 
value of some observable O in this system is given by 

{0) = - J V<p 0[0]exp(-%]) , where Z = J exp(-S[<p]) (1) 

is the partition function, and Dip is the sum over all closed paths <p(t) G S 1 in some period T, 
i.e. t £ [0, T] and (p(Q) = <p(T). The set of all these paths is naturally divided into disjoint 
subsets, which are characterised by the winding number 



Q = ±- [ T dt<p^% , (2) 
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which represents in this case the topological charge. Continuous path deformations cannot change 
Q, hence these subsets are indeed topological sectors. 



Topological sectors also occur in a variety of gauge theories [T]. Let us consider gauge 
configurations in a Euclidean space with periodic boundary conditions (a torus). If they are 
split into topological sectors, the characteristic topological charge Q is also denoted as the 
Pontryagin index. Two examples are 



where F pu is the field strength tensor, and F pv := e pup<J F pa . Gauge configurations can be con- 
tinuously deformed only within a fixed topological sector, hence the functional integral splits 
into separate integrals for each Q G Z. 

Let us now address such a gauge theory in the presence of chiral fermions, i.e. massless 
fermions with a Dirac operator D that anti-commutes with 75, D75 + 75Z? = 0. In this case the 
zero modes of the Dirac operator have a definite chirality ±1. 

For such a Dirac operator, in a given gauge background, we denote the number of zero modes 
with chirality +1 (—1) as n + (n_). Their difference is the fermion index 



The Atiyah-Singer Index Theorem \2\ states that for any gauge configuration, this index coincides 
with the topological charge 



2. Lattice regularisation 

The lattice regularisation reduces the (Euclidean) space to discrete sites x, which are 
separated by some finite lattice spacing a. The latter implies an UV regularisation of the 
corresponding quantum field theory. Matter field variables are now defined on each lattice site, 
e.g. ty x , ^> x for fermion fields, while gauge fields can be formulated as compact link variables 
U x ,fi £ { gauge group } . It is a great virtue that this formulation is gauge invariant even on the 
regularised level, so in this approach no gauge fixing is needed. 

A priori there are no topological sectors anymore in the lattice regularised system; all 
configurations can now be continuously deformed into each other. Still, the desired connection to 
the continuum theory motivates the attempt to introduce somehow (the analogue of) topological 
sectors also on the lattice. A number of suggestions appeared in the literature, often with a 
somewhat questionable conceptual basis. A clean formulation emerged only at the very end of 
the last century, based on chiral lattice fermions. Their lattice Dirac operator D cannot simply 
anti-commute with 75 due to the notorious doubling problem of lattice fermions [3] , but it may 
obey the Ginsparg-Wilson Relation (GWR), which reads (in its simplest form) 




(3) 



v ;= n_ — n + € "Z . 



(4) 



Q = v . 



(5) 



-D75 + = aD^ 5 D . 



(6) 



This still guarantees a lattice deformed — but exact — version of the chiral symmetry [1]. The 
latter also implies that the corresponding lattice Dirac operator has exact zero modes with a 



definite chirality, as in the continuum. Hence we can adopt the Index Theorem [5] and define 
the topological charge of a lattice gauge configuration as Q := v. 

We remark that random lattice gauge configurations always occur with n+ = or ra_ = 0; 
configurations with a cancellation in the lattice fermion index also exist (the free fermion is an 
example), but their probability measure seems to vanish 

3. Monte Carlo simulation 

Observables in quantum gauge theory can be evaluated beyond perturbation theory by means 
of Monte Carlo simulations. The idea is to use a sizeable set of gauge configurations [U] 
(consisting of link variables all over the lattice volume), which are generated randomly with 
the probability distribution 

p[U] = det D[U] exp(-5 gaugc [C/]) . (7) 

Here we assume a fermion action, which is bilinear in the Grassmann valued spinor fields ^, ^. 
Their functional integration T>fyT>*$> is carried out already, giving rise to the fermion determinant 
det D[U]. 

The summation over this set of configurations yields a numerical measurement of expectation 
values (...), in particular of n-point functions. These results obviously come with some statis- 
tical error (since the available set of configurations is finite), and a systematic error (e.g. due 
to the finite lattice spacing a, which usually requires a continuum extrapolation). Both can be 
estimated and reduced if necessary by extended simulations. On the other hand, we stress again 
that the result is fully non-perturbative; we deal with the complete action in the exponent, i.e. 
we capture directly the given model at finite interaction strength. 

Practical algorithms for the generation of gauge configurations (with the given probability 
distribution) perform local updates, i.e. in one step a configuration is modified just locally. 
Iterating such steps many times leads to a (quasi-)independent new configuration, to be used 
for the next measurement. Changing a gauge configuration drastically in a single step is also 
conceivable in principle, but in practice such algorithms tend to be inefficient 

A problem with a sequence of local updates is, however, that it hardly ever changes the 
topological sector — although one should do so frequently in order to sample correctly the 
entire space of configurations. This problem is particularly striking in the attempts to simulate 
QCD with chiral quarks; the JLQCD Collaboration performed very extensive 2-flavour QCD 
simulations of this kind [6] — which led to interesting results — but the Monte Carlo histories 
were always confined to the trivial topological sector of charge Q = 0. 

Most QCD simulations with dynamical quarks involve a non-chiral lattice quark formulation, 
since Ginsparg- Wilson fermions are tedious to simulate. In particular Wilson fermions (and 
variants thereof) have the disadvantage of additive mass renormalisation, but the problem with 
sampling the topological sectors is less severe so far. However, that property depends on the 
lattice spacing; typical values that have been used in the past are a ~ 0.05 . . . 0.1 fm. Once one 
tries to proceed to even finer lattices, the problem of confinement of the Monte Carlo history to 
a single topological sector is expected to show up also in this formulation [7jH 

The same holds for the specific lattice field configurations which are exactly on a topological boundary, so we 
can ignore them. 

2 Cluster algorithms are a counter example for certain spin models, but no efficient application to lattice gauge 
theories is known so far. 

3 To be more explicit: any algorithm has to obey "detailed balance", i.e. the transition probabilities of some 
configuration Ci to C2 and vice versa have to match the probability ratio for these configurations to occur 



So we have to address the question how to handle Monte Carlo simulations if the history 
tends to be trapped for a very long (computing) time, i.e. for many, many update steps, in one 
topological sector. What are then the prospects for measuring some n-point function, or the 
topological susceptibility 

Xt •= \ ((Q 2 ) - (Q) 2 ) {V = volume) , (8) 

which actually require the summation over a variety of topological sectors, with suitable 
statistical weights? 

This is a delicate and highly relevant issue. Here we address it in a toy model study of the 
2-flavour Schwinger model, which we simulated [5] with dynamical overlap hypercube fermions; 
this is one version of chiral lattice fermions [S], with a Dirac operator that solves the GWR JU||j 
We designed and applied a variant of the Hybrid Monte Carlo algorithm, which is particularly 
suitable for this type of lattice fermions [5] . 



4. The Schwinger model 



The Schwinger model [TT] represents Quantum Electrodynamics on a plane (QED2). It is a 
popular toy model; in particular it shares with QCD the property of fermion confinement |12j 
(although the gauge group is Abelian) and the presence of topological sectors, see eq. ©. On 
the other hand there are qualitative differences, such as the absence of a running gauge coupling 
in the Schwinger model. In the continuum its Lagrangian can be written as 

Ffiv {x 

)• (9) 

We are interested in the case of Nf = 2 degenerate fermion flavours of mass m<ft where Ref. 
|13| made the following predictions: 

chiral condensate £ := = 0.388 . . . (mg 2 ) 1/s , (10) 

pion mass M n = 2.008 . . . (m 2 g) 1/3 . (11) 

As in 2-flavour QCD, a "meson" singlet and a triplet emerge, the former (latter) being massive 
(massless) in the chiral limit m — > 0, cf. eq. (|11|) . Referring to this analogy, and in agreement 
with the literature, we denote the triplet as "pions" . Its emergence in 2 dimensions might appear 
somewhat surprising; the theoretical background of these "quasi-Nambu-Goldstone bosons" was 
first discussed in Ref. 



5. Numerical measurement at fixed topology 



We simulated the 2-flavour Schwinger model at/3:=l/g 2 = 5[8]. This implies smooth gauge 
configurations (mean plaquette value ~ 0.9). Also the "meson" dispersion relations confirm 
that lattice artifacts are tiny [8], hence we can confront our results directly with the continuum 
predictions (|10p . (jlip . without really needing a continuum extrapolation. On the other hand, 
finite size effects are significant, and they are in fact necessary for our discussion of topology 
dependent observables. We simulated on L x L lattices of sizes L = 16, 20, 24, 28, 32 with 

(Boltzmann weights), p(Ci — > C2)/p(C2 — > Ci) = exp(S[Ci] — S[C2])- The boundaries between topological 
sectors are surrounded by zones of high action, i.e. low probability. As the lattice spacing a is reduced, their 
weight p(C) decreases with a high power of a [7]. Hence a sequence of small update steps will rarely tunnel 
through such a boundary. 

4 In this formulation we insert an improved kernel into the overlap formula, instead of the Wilson kernel of the 
standard overlap operator [10]. The virtues include an improved locality and scaling behaviour, and approximate 
rotation symmetry [9]. 
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Table 1. Statistics for lattice size L = 16 (above) and L > 16 (below) and various fermion 
masses m, in distinct topological sectors (y is the fermion index). Multiple starts of Monte Carlo 
histories were necessary to get access to various topological sectors, since topological transitions 
were very rare. 

fermion masses in the range m = 0.01 . . . 0.24 (in lattice units). Our statistics is displayed in 
Table E 

Let us first address the Dime spectrum. All the eigenvalues of a lattice Dirac operator (before 
adding the mass), which obeys the GWR ©, are located on the circle in the complex plane 
with centre and radius 1/a, as illustrated in Fig. [TJ This confirms that the zero modes are exact, 
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Figure 1. The spectrum of a lattice Dirac operator that fulfils the GWR © is located on a 
circle in the complex plane. In the continuum limit a — > it turns into the imaginary axis. 

and we have mentioned before that their fermion index is identified with the topological charge, 
Q = v = n_ — n + 

In this study we could evaluate the complete Dirac spectrum for our lattice configurations 
5 Actually throughout this study only the absolute value \v\ matters. 



(which is not feasible in 4 dimensions, except for tiny lattices). To make this spectrum compatible 
with the continuum formulation, we map it stereographically onto the imaginary axis [15]. Based 
on the eigenvalues Aj that we obtain after this mapping, we obtain the chiral condensate 




The sum can be computed for each configuration, but expectation values can only be measured 
within fixed topological sectors. Table [1] shows that topological transitions are indeed so rare 
that the entire space of configurations is not well sampled, but specific sectors are explored well. 
Hence we measure results for the expectation values of the chiral condensate at specific values 
of \u\, 




In the last expression we split off the zero mode contribution to Sij/i, which dominates at small 
mass m (and v ^ 0), and we denote the rest as E\ v \. Numerical results are shown in Fig. 




Figure 2. Numerical results for the chiral condensate at [topological charge) = \v\, after 
subtracting the zero mode contribution, cf. eq. ([T5j) . at L = 16 (on the left) and at m = 0.01 
(on the right). 

It is a generic property of stochastic Hermitian matrices (such as 75-D) that zero eigenvalues 
repel the low-lying non-zero modes. This suggests the inequality 

£0 > ei > £2 • • • (14) 

at fixed m and V, which is confirmed consistently by the plot in Fig. [2] on the left. Moreover 
the plot on the right shows that 

> £i(V 2 ) for Vx > V 2 , (15) 

which is less obvious: in a larger volume more eigenvalues cluster near zero, which supersedes 
the pre-f actor 1/V. 

The rest of this article is devoted to tests of three different methods for approximately 
extracting "physical" quantities (i.e. quantities which are properly summed over all topological 
sectors), based on measurements in a few specific sectors. 



6. Gaussian evaluation of the topological susceptibility 



We first assume a Gaussian distribution of the topological charges - - this is certainly 
reasonable, for instance precision tests in SU(3) pure gauge theory revealed at most tiny 
deviations from this behaviour |16j . It implies that the chiral condensate is composed as 

Parity symmetry assures that (u) = 0, hence the topological susceptibility simplifies to 

X, - ^ . (IV) 

In most volumes we have data for Eo . . . £q, i.e. up to some maximal topological charge Q. 
Thanks to inequality (jl4|) all the higher charge contributions — for \u\ > Q — are bounded as 

^-L <£ M < ^-L + eg . (18) 
mv 1 1 mv 

Hence for a given value of the susceptibility xt the sum in eq. (fl~6j) can be performed, up to a 
uncertainty which affects £ only mildly, since Ei„i for high charges contribute only little. 

In two volumes, L = 24 and 28, some Ei„i data are missing for \v\ < Q (see Table[T]); in these 
cases we can again fix a minimal and a maximal value for e\ v \ , this time based on inequality (|15p 
and the results in the next smaller and next larger volume. 

So we can probe any ansatz for \t and compute the corresponding value of £ up to a modest 
uncertainty. We require the result to agree (within errors) with the prediction (|10p . In this way 
we determine xt- Fig- El shows the results for L = 16 and m = 0.01, 0.03, 0.06 — for higher 
masses the assumption m -C g — 0.45, which is needed for the prediction (|10p . seems to fail. 
Since the theory refers to infinite volume, we expect the result to improve for increasing m (i.e. 
for shorter correlation length) within the allowed range. 




Figure 3. Our results for the topological susceptibility xt based on Gaussian summation. On 
the left: in a fixed volume V = 16 2 our data are compatible with an interpolation Xt( m ) \f™i 
and for increasing mass (i.e. reduced finite size effects) we approach the conjecture (|19p . On 
the right: at fixed m = 0.01 we obtain results for xt, which seem to follow a behaviour linear in 
1/V. 



The result is compared to a QCD-inspired conjectured of Ref. [T7] (for Nf flavours, in a large 
volume) , 



1 

Xt 



m 



E(N f = l,m = 0) + 



Xt(N f = 0) 



(19) 



The first ingredient has been computed analytically, E(iV/ = l,m = 0) ~ 0.16 <7 |14j . and the 
quenched susceptibility Xt(Nf = 0) ~ 0.000332 has been measured numerically [T5]. Fig. 
confirms that the corresponding curve approaches the fit through our values for increasing m. 

Alternatively we fix the mass m = 0.01 and compare the results for L = 16 ... 32 (plot in 
Fig. on the right). In our largest volume, V = 32 2 , we obtained xt = 0.000341(4), which is 
close to the value given by conjecture (fTUj) . \t = 0.000332. An infinite volume extrapolation of 
our data, however, leads to a smaller susceptibility of xt = 0.000226(5). 

7. Correlation of the topological charge density 



A drawback of the method in Section 6 is that a known reference quantity is needed (here 
it was £), and results in various topological sectors are required. This is not the case for an 
approach suggested in Ref. [19] . which derived a "model independent formula" for the correlation 
of the topological charge density pt in one sector, 



lim 



(p t (x)pt(0))|, 



1 V 2 

yXt+y5+0{V- 



(20) 



For tests in 2-flavour QCD we refer to Ref. [2D]. (The original formula even includes a correction 
for a possible deviation from a Gaussian distribution of the topological charges, which we 
neglect.) In order to justify the assumptions in the derivation of this formula, we have to 
assume a large expectation value (v 2 ) = Vxt, and a small ratio \u\/(i> 2 ). 

As an example, we show in Fig. [5] the corresponding correlation at L = 16, v = and 
various masses. Numerically the density was computed from the simplest lattice version of 
Pt = £12-^12 (this is not problematic in the current setting, where we are always dealing with 
smooth configurations) . 
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Figure 4. The correlation of the topological charge density for L = 16, v = and 
m = 0.01, 0.03, 0.06. The statistical noise does not allow us to resolve a plateau value at 
large distances — overcoming this problem would require a huge statistics. 

At large distances one should find a plateau value, which would then yield xt- In particular 
for m = 0.06, where the maximal distance might be sufficient to see the asymptotic behaviour, 



we expect (based on the data and the conjectured formula in Section 6) a plateau value of 
Xt/V ~ —6 • 1CP 6 . However, our statistical errors are of 0(1CP 5 ), so in order to clearly resolve 
this plateau we would need about 50 000 to 100000 configurations (cf. Table [I]). We conclude 
that the applicability of this method requires unfortunately a very large statistics. 



8. Approximate topological summation of observables 



We now proceed to the main approach in this study. It is a method that does not require a 
known input observable either (as in Section 7), but measurements in various topological sectors 
and volumes are needed. In fact this is the input which is usually accessible. Then one tries to 
extract a (topologically summed) observable (O) by employing the approximation formula 

<0) M «<0> + ^1-^ (c = const.). (21) 

This formula has been derived first for the pion mass in QCD |21| . but it applies generally 
to observables in a field theory with topology [8]. As in Section 7 one assumes a Gaussian 
distribution of the topological charges, and a large value of Vxt, as well as a small ratio 
\u\/(u 2 }, are favourable for the validity of the approximations involved in the derivation. This 
approximation formula could be truly powerful in QCD and elsewhere, but it has never been 
tested before. 



8.1. Application to the chiral condensate 



Let us apply formula (|2ip to the chiral condensate. It is convenient to modify the notation, 

E^E-^ + f 2 ^, A = --, 5 = -4- (22) 
V V 2 xt Xt 

The unknown quantities are £, A and B, and we are ultimately interested in £ and xt = A/B. 
They can be determined (in the framework of this approximation) by numerical results for some 
E„: 

• At fixed m and V, we can determine B, for instance from Eo and Ex- 

• If we keep m fixed but consider two volumes, V\ ^ we can further determine A, e.g. 
based on Eo- 

In total, it takes (at least) three E y values, involving two volumes, to obtain results for E and 
Xt- 

We follow this sequence of steps and start with the determination of B. If we use as our 
input the measurements in the topological sectors with \u\ = k, I (at fixed m and V), we denote 
the result as Bk/, 

1 Efc — Yj£ 1 £k — Eg 

v Bk * = v l^¥ = MkTT) + v W^p ■ (23) 

The semi-classical term, l/(m(k + £)), tends to vary strongly for different choices of k and I. 
Ideally the quantum effects should render the results for B^/ similar again. As an example, we 
show in Fig. [5] results for B\. o at L = 16, m = 24. In fact the non-perturbative results are much 
more stable in k than the semi-classical contributions alone. Hence the first consistency test is 
passed well. 



We proceed to the determination of A, and therefore of E, based on Eo measurements in two 
volumes with sizes (Li,L2). Here we consider m = 0.01 and we give two examples: 
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Figure 5. Results for the auxiliary variables Bk : o (cf. eqs. ([22]) . ([23]) ) at L = 16, m = 0.24, 
for k = 1,...,4. We see that the numerical results are quite consistent, in contrast to the 
semi-classical values. 



• The S values in (Li,L 2 ) = (16,32) yield E = 0.0199(7). 

• The E values in (Li,L 2 ) = (20,32) yield E = 0.0207(12). 

Thus the consistency looks fine again, but these results are far below the prediction ()10p . 
E = 0.04888 (in an infinite volume). In this case, our results are strongly affected by finite 
size effects, which is not surprising: for the given fermion mass, the correlation length in infinite 
volume (given by eq. (fTT]) ) would be £ = l/M 7T y =OQ ~ 14. The relatively small boxes enhance 
the Dirac eigenvalues |Aj|, such that E decreases. 

So the mass m = 0.06 should be more promising, where theory predicts £ ~ 4.3. Here we only 
have data in [L\,Li2) = (16,32), for \v\ = 0, 1, so we cannot repeat the above consistency tests. 
Nevertheless we can evaluate xt = 0.00118(30) (which is just compatible with the conjecture 
(|19p . xt — 0.00146). We further insert our most reliable result for B, namely B\^ measured in 
L = 32, and arrive at a result for E, which is indeed close to the theoretical prediction (fTT]) . 

Enumerical = 0.0883(69) , S thcory = 0.0888 . (24) 
8.2. Application to the pion mass 

Let us also test the approximate summation formula (|2ip by applying it to the pion mass. 
As we mentioned before, this was the original idea of Ref. |21| (though that work referred to 
QCD). We re- write approximation ([2T]) in the notation analogous to ([22]) . 

M ffiM « M n - ^ + ^u 2 . (25) 

However, we now adopt a strategy which differs from the previous consideration of E: at fixed 
m we determine the three unknown parameters A, B, M n directly by a least-square fit for some 
set of numerical M^\ v \ values. 

• For m = 0.01 we have in total 11 measurements of (see Table [1]), and we include 
the most promising ones. We need at least two volumes, so we take the largest two with 



(L\,L2) = (28,32). Moreover we only include the topological sectors with \u\ < 1, which are 
favourable for the condition that \v\/(u 2 } should be small. This leads to 

M w ,q M w ,i ) fit 
L = 28 : 0.146(4) \ ^ M n = 0.073(25) , (26) 

L = 32 : 0.05(1) 0.160(8) J 

which matches well the theoretical prediction, M n = 0.071 (albeit with a large error). 

• We proceed to m = 0.06, where we only have data for (Li, L2) = (16, 32). Hence we have 
less choice in this case, but the finite size effects are less severe. Again we include the results for 
\u\ < 1, which corresponds to four input measurements this time, and we arrive at 

ik^o M ff ,i ) fit 
L = 16 : 0.041(1) 0.271(4) > M n = 0.233(8) . (27) 

L = 32 : 0.23(1) 0.232(7) J 

Also this result agrees well with the theoretical pion mass, M n = 0.235, and this time also the 
uncertainty is modest. 

9. Conclusions 

We have addressed a quite generic problem of lattice simulations in gauge theories with 
dynamical ( quasi- )chiral fermions. The Monte Carlo histories of such simulations tend to get 
trapped in one topological sector for a very long (simulation) time, i.e. over many update steps 
of the lattice gauge configuration. A conceptual issue that one has to address in this situation 
is ergodicity, a property which is compulsory for a correct algorithm. Here we studied a more 
practical question: how can we evaluate the expectation value of some observable (O), when 
only numerical measurements restricted to a few topological sectors, {0) u , are available? 

The dominant subject in contemporary lattice simulations is QCD with dynamical quarks. 
Here the problem of topological restriction is most striking when one deals with chiral lattice 
quarks (of overlap [10] or Domain Wall [22] type) , which solve the Ginsparg- Wilson Relation (eq. 
([6]) or generalisations thereof). The use of Wilson type quarks is more widespread because they 
are much faster to simulate, though plagued by additive mass renormalisation and problems 
related to operator mixinglflQ 

Here the aforementioned topological problem is less severe so far, but it is expected to show 
up as well when simulations will be carried out on finer and finer lattices, say with lattice spacing 
a < 0.05 fm. This renders the lattice QCD formulation more and more continuum-like, which 
is in general welcome, but it also makes it more difficult to change the topological sector. 

This problem is not manifest in a very large volume, where {0) v is the same for all indices 
v (this property agrees with approximations (|20|) . (|21|) ). However, to suppress the topological 
dependence and other finite size effects, the volume has to be large compared to the correlation 

6 Also that problem is avoided by the use of Ginsparg- Wilson fermions [23] . 

7 For completeness we add that "staggered fermions" are widespread as well in lattice QCD. They are also quick 
to simulate, and they do not suffer from additive mass renormalisation, but the number of flavours is not flexible. 
Therefore it is now popular to take the fourth root of the fermion determinant (cf. eq. (0), which formally 
corresponds to a single flavour, but this is harmful for locality, which is conceptually important. The question if 
this is a reason to worry in practice is highly controversial. In any case, neither Wilson nor staggered fermions 
do provide a sound definition of the topological charge since there is no well-defined fermion index, in contrast to 
Ginsparg- Wilson fermions [S]. Hence one has to refer to some rather hand- waving definition in these cases. 



length, which is given by the inverse pion mass, £qcd ~ 1-4 fm <C La. But when a is very small, 
this requires a huge lattice size L, which makes simulations again very tedious. 

As a way out, the use of open boundary conditions in the Euclidean time direction has re- 
cently been advocated, so that topological charge can gradually flow in or out of the volume 
during a simulation |24j . In our study, however, we stay with periodic boundary conditions 
for the gauge fields, which guarantee that the topological charge is always integer, along with 
(discrete) translation invariance. As a toy model we considered the Schwinger model with two 
light, degenerate flavours, which were represented on the lattice by dynamical overlap hypercube 
fermions. In a set of small or moderate volumes, this only enabled measurements inside some 
specific topological sectors. In order to establish a link to the "physical" quantities, we tested 
three methods to approximate the topological summation: 

• The confrontation of a Gaussian summation with a known observable allows us to fix the 
topological susceptibility xt- This method is robust, but it requires a known input quantity. 
This is available in the 2-flavour Schwinger model [13] (we used the chiral condensate), but 
not in general. 

• Next we tested a method to evaluate xt based on the correlation function of the topological 
charge density [19] . More precisely, one searches for an asymptotic plateau of this correla- 
tion at large distances, which should amount to — Xt/V ( a t v = 0). Unfortunately this value 
tends to be tiny for realistic settings, hence its resolution requires a very large statistics. 

• Our main goal was the test of an approximate summation formula given in Ref. [21], 
which could provide a "physical" result (O), using only measurements of some topologically 
restricted observables (0) u as an input — for various values of in at least two volumes. 
This method is potentially powerful, but it has never been tested before. 

Our results suggest that it may work, if the assumptions used in the derivation of this 
formula are reasonably well justified. In particular, Vxt = (^ 2 ) should be "large", but it is 
difficult to predict explicitly what this means. In our settings this quantity was always below 
0.5, but nevertheless we found decent (though not very precise) results for the topologically 
summed chiral condensate and pion mass. This observation is encouraging for applications 
in QCD simulations with dynamical quarks. 
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